A synthetic electric force acting on neutral atoms 
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Electromagnetism is a simple example of a gauge theory where the underlying potentials - the 
vector and scalar potentials - are defined only up to a gauge choice. The vector potential generates 
magnetic fields through its spatial variation and electric fields through its time-dependence. We 
experimentally produce a synthetic gauge field that emerges only at low energy in a rubidium 
Bose-Einstein condensate: the neutral atoms behave as charged particles do in the presence of a 
homogeneous effective vector potential. We have generated a synthetic electric field through the 
time dependence of an effective vector potential, a physical consequence even though the vector 
potential is spatially uniform. 



Gauge theories play a central role in modern quantum 
physics. In some cases, they can be viewed as emerging 
as the low energy description of a more complete the- 
ory [TJ [2]. Electromagnetism is the best known gauge 
theory and its gauge fields are the ordinary scalar and 
vector potentials. Magnetic fields arise only from spa- 
tial variations (curl) of the vector potential, while elec- 
tric fields arise from both time variations of the vector 
potential and gradients of the scalar potential; these po- 
tentials are defined only to within a gauge choice. For 
a charged particle in such potentials the canonical mo- 
mentum (the variable canonically conjugate to position) 
and the mechanical momentum (mass times the veloc- 
ity) are not equal. Our experiments [3] have realized a 
particular version [4] of a class of proposals [5HTD] to 
generate effective vector potentials for neutral atoms via 
interactions with laser light, and have created synthetic 
magnetic fields [11] on the success of simulating charged 
condensed matter systems with neutral atoms [12] [13]. 
Here we demonstrate the complementary phenomenon: 
the generation of a synthetic electric field from a time- 
dependent effective vector potential. Additionally, we 
make independent measurements of both the mechanical 
momentum and canonical momentum, where the latter 
is usually not possible. 

Introduction The electromagnetic vector potential A 
for a charged particle appears in the Hamitonian H = 
(Pcan — #A) 2 /2ra, where p can is the canonical momen- 
tum, q is the charge, and m is the mass. (The quan- 
tity p can — qA — mv is the mechanical momentum for 
a particle moving with velocity v.) We recently demon- 
strated a technique to engineer Hamiltonians of this form 
for ultra-cold atoms, and prepared a Bose-Einstein con- 
densate (BEC) at rest with an effective vector potential 
A — A x x constant in time and space [3 , correspond- 
ing to E = B = 0. In Ref. [11 , we made A depend 
on position giving B = Vx A / but E = 0. Here 
we add time dependence to a spatially uniform vector 
potential A(t) = A(t)x, generating a synthetic electric 
field E(t)x = —dA/dt. The resulting force is distinct 
from that due to fields arising from gradients of scalar 



potentials 0(r) like gravity or an external trapping po- 
tential. A revealing analog to which we will return is 
that of an infinite solenoid of radius ro as pictured in 
Fig. 1A: a magnetic field B = Bz exists only inside the 
coil, however, a non-zero cylindrically symmetric vector 
potential A = Br^j2r in the symmetric gauge extends 
outside the coil. Far from the coil, A is nearly uniform, 
analogous to our uniform effective vector potential. 

We experimentally synthesize electromagnetic fields 
for neutral atoms by illuminating an atomic 87 Rb BEC 
with two intersecting laser beams (Fig. 2A) that cou- 
ple together three atomic spin states within the 5S elec- 
tronic ground state. The three new energy eigenstates, 
or "dressed states," are superpositions of the uncou- 
pled spin and linear-momentum states and have energy- 
momentum dispersion relations that can be dramatically 
different from those of uncoupled atoms. The atoms in 
a dressed state act as particles with a single well-defined 
velocity v, which is the population- weighted average of 
all three spin components. 

The dispersion relation of the lowest energy dressed 
state changes near its minimum, from p 2 /2m to (p — 
Pmin) 2 /2m* (as shown in Fig. 2C), where the minimum 
location p min plays the role of qA. In addition, the 
mass m is modified to an effective mass m*, m* > m, 
and both p m i n and m* are under experimental control 
(although not completely independently). We identify 
Pmin = #*A*, the product of an effective charge q* and 
an effective vector potential A* for the dressed neutral 
atoms. As we change A* in time, we induce a synthetic 
electric field E* = -dA*/dt, and the dressed BEC re- 
sponds as d{m*v)/dt = — V0(r) + <?*E*, where v is the 
velocity of the dressed atoms and m*v = p can — g*A*. 
Here, A(m*v) = — <?*(Aj — A*) is the momentum im- 
pulse imparted by the synthetic electric force g*E*. 

We study the physical consequences of sudden tempo- 
ral changes of the effective vector potential for the dressed 
BEC. These changes are always adiabatic such that the 
BEC remains in the same dressed state. We measure 
the resulting change of the BEC's momentum, which is 
in complete quantitative agreement with our calculations 
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FIG. 1. Schematic of the electric field generated by a time varying vector potential. (A) Emulated system, showing the electric 
current flowing counter clockwise in the infinite solenoid (black coil) with radius ro and the real magnetic field B only inside 
the solenoid. The blue lines represent the vector potential A. A charged particle (red dot) located far from the coil experiences 
a nearly uniform A. (B) Calculated time response of the synthetic vector potential and electric field for neutral atoms in our 
first measurement (see Fig. 3). The calculation includes the known inductive response time of the bias field Bo, which sets the 
detuning, and the calibration of detuning to vector potential shown in Fig. 2D. 



and constitutes the first observation of synthetic electric 
fields for neutral atoms. 



Experimental realization Our system, depicted in 
Fig. 2A, consists of a F = 1 87 Rb BEC with about 
1.4 x 10 5 atoms initially at rest [HJIT5; a small physical 
magnetic field Bo Zeeman-shifts each of the three mag- 
netic sublevels mp = 0, ±1 by an amount £^o,±i 7 but pro- 
duces no Lorentz force on our charge-neutral BEC. Our 
experiment is performed in the linear Zeeman regime, 
i.e., E-i & —E+i & g\±BB§ ^> \Eq\, with linear and 
quadratic Zeeman shifts huz = (E-i — E+i)/2 « h x 
2.32 MHz and -he = Eq-{E_ 1 +E +1 )/2 « -ft x 784 Hz. 
In addition, a pair of A = 801 nm laser beams, inter- 
secting at 90° at the BEC, couples the mp states with 
strength Qr. These Raman lasers differ in frequency by 
Aul ~ ooz and we define the detuning from Raman reso- 
nance as 5 = Aujl—ojz- In our experiments hQf> ~ lOEp 
and \HS\< WE L , where E L = H 2 k 2 L /2m and k L = V^tt/X 
are the appropriate units of energy and momentum. 

When the atoms are rapidly moving or the Raman 
lasers are far from resonance (fc^v or S ^> £Ir), the cou- 
pling lasers hardly affect the atoms. However, nearly 
resonant, slowly moving atoms are strongly dressed: 
the three uncoupled states transform into three new 
dressed states. The spin and linear- momentum state 
\k x ,rriF = 0) is coupled to states \k x — 2ki Jl mp — +1) 
and \k x + 2/c^,^f = —1), where hk x is the momentum 
of | rap = 0) along x, and 2hki J x is the momentum dif- 
ference between two Raman beams. In the frame rotat- 
ing at Aujl under the rotating wave approximation, the 
Hamiltonian matrix H(k x )/h for motion along x is 



/ ^(k x + 2k L ) 2 -d n R /2 \ 

n R /2 ^kl-e n R /2 

V Q R /2 J^(k x -2k L ) 2 +8 J 

Diagonalization for each k x yields three dressed eigen- 
states, each of which is a superposition of the uncoupled 
states [3], with energies Ej(k x ) shown in Fig. 2C (gray 
for uncoupled states, colored for dressed states); we focus 
on atoms in the lowest energy energy-momentum (red) 
curve. When the atoms' energy (interaction and kinetic) 
is small compared to the « IOEl energy difference be- 
tween the curves-as it is in our experiments-the contact- 
interacting atoms will remain within the lowest energy 
dressed state [4]. Such low energy particles do not re- 
veal their spin and momentum components. In the low 
energy limit, dressed atoms have a new effective Hamil- 
tonian for motion along x , H x = (hk x — q*A%) 2 /2rrf 
(motion along y and z is unaffected), where we have cho- 
sen a gauge such that the momentum of the mp = 
component hk x = p can is the canonical momentum of the 
dressed state. Our approach introduces the term 
which depends on the detuning 5 as shown in Fig. 2D, 
and an effective mass m*. The red curve in Fig. 2C shows 
H x for <?*A* > and the arrow indicates the canoni- 
cal momentum of a BEC in the lowest energy dressed 
state, p can = q*A* x . Although this dressed BEC is at rest 
(v = dH x /dhk x = 0, zero group velocity), it is composed 
of three bare spin states each with a different momentum, 
among which the momentum of |rap = 0) is hk x = p C an- 
Although none of its three bare spin components has zero 
momentum, the BEC's averaged momentum is zero. 

We measure p can by first removing the coupling fields 
and trapping potentials and then allowing the atoms to 
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FIG. 2. Experiment setup for synthetic electric fields. (A) Physical implementation indicating the two Raman laser beams 
incident on the BEC (red arrows) and the physical bias magnetic field Bo (black arrow) . The blue arrow indicates the direction 
of the synthetic electric field E* . (B) The three wif levels of the F — 1 ground state manifold are shown as coupled by the 
Raman beams. (C) The dressed state eigenenergies as a function of canonical momentum for the realized coupling strength of 
MIr = 10.5El at a representative detuning HS — —IEl (colored curves). The grey curves show the energies of the uncoupled 
states, and the red curve depicts the lowest energy dressed state in which we load the BEC. The black arrow indicates the 
dressed BEC's canonical momentum p ca n = q* A* , where A* is the vector potential. (D) Measured vector potential from the 
canonical momentum. 



freely expand for a t = 20.1 ms time-of- flight (TOF). 
Because the three components of the dressed state 
{\k x ,rriF = 0),|/c x ^2&L,raF = ±1)} differ in momen- 
tum by ±^2/cl, they quickly separate. 

Using a combination of adiabatically applied Raman 
and rf (radio-frequency) fields we transfer (load) the BEC 
initially in \rriF = —1) into the lowest energy dressed 
state with a A* = A*x (see Ref. [3] for a complete tech- 
nical discussion of loading), i.e., with a canonical momen- 
tum p can = q*A*. We measure as the momentum of 
\™>f = 0); Fig. 2D shows how the measured and predicted 
A* depend on the detuning S. With this calibration, we 
use S to control the effective vector potential. 

Synthetic electric field Electric fields result either 
from spatially varying scalar potentials or from time 
varying vector potentials. An example of the latter case 
is illustrated in Fig. 1A for an infinite solenoid. The 
magnetic field B exits only inside the coil; when the cur- 



rent is constant the electric field E is everywhere zero 
but the symmetric gauge vector potential is zero only 
at the solenoid center. When the current is changed 
in a time interval At the magnetic field inside the coil 
changes to a final value B(At) , and the exterior vec- 
tor potential changes by AA = [B(At) - B(0)] rg</>/2r. 
Here, a charged particle on the y axis feels an electric 
field —(dA/dt)x during the time At, leading to a me- 
chanical momentum kick Ap = — qAAx, even outside 
the solenoid. 

We realize a synthetic electric field E* for ultracold 
neutral atoms by changing the effective vector potential 
from an initial value A* to a final value Aj. We perform 
two types of measurements of E*: (1) we remove the 
trapping potential right before we change A* such that 
E* can accelerate the atoms unimpeded, and measure the 
resulting change from zero of the BEC's velocity; and (2) 
after changing A*, we allow the trapped, dressed atoms 
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FIG. 3. Change in momentum from the synthetic electric field 
kick. This figure displays three distinct sets of data shown 
by red and blue circles and black crosses. The data were 
obtained by applying a synthetic electric field by changing 
the vector potential from q*A* (between +2/i/cl and — 2/i/cl) 
to q*A*f. Circles indicate data where the external trap was 
removed right before the change in A* : for red (blue) symbols 
q*A*f = —2hkj J (-\-2hki l ). The black crosses, more visible in the 
inset, show the amplitude of canonical momentum oscillations 
when the trapping potential was left on after the field-kick. 
The typical uncertainties are also visible in the inset. The 
grey line is a linear fit to the data (circles) yielding slope 
—0.996 =b 0.008, where the expected slope is —1. 



to oscillate in their harmonic potential, and measure the 
magnitude and frequency of the oscillations of both the 
canonical momentum and velocity. 

We prepare our BEC at rest with an initial vector po- 
tential A*x. In the first measurement, we change the 
detuning S in 0.8 ms, leading to a change of the vec- 
tor potential to A^ after which the Raman coupling is 
turned off in 0.2 ms. As A* changes from A* to A*f, 
a synthetic electric field E* appears and accelerates the 
BEC. We image the atoms after a TOF and record their 
final velocities. Figure 3 shows the momentum Ap im- 
parted to the atoms by E* as a function of the vector 
potential change q*(A*j— A*), denoted by red (blue) sym- 
bols with q*A* f /h = 2fc L (— 2fc L ). We have chosen a suf- 
ficiently large final detuning HS = ^60El such that the 
final atomic state is a nearly pure spin state, \mp = +1) 
for q*A* f /H = 2k L or \m F = -1) for q*A* f /h = -2k L . 
For this essentially undressed final state, m* = m and 
the imparted momentum is Ap = m*v = mv, equal to 
the change in mechanical momentum. We performed a 
linear fit Ap = Cq*(A*j — A*) to the data and obtain 
C — —0.996(8), in good agreement with the expected 
Ap=-q*(A* f -A*). 

In the measurement described above, the final value 
of A* = A*j resulted from a detuning so large that the 
final-state atoms were essentially undressed. In a second 
measurement, we examined the time evolution of atoms 
that remain trapped and strongly dressed. Here we again 



changed A* from A* to A J in about At = 0.3 ms but left 
the dressed BEC in the harmonic confining potential for a 
variable time before abruptly removing the confining po- 
tential and the Raman fields for a TOF. As the BEC os- 
cillated in the trap, we monitored the out-of-equilibrium 
canonical momentum p ca n- It is our access to the inter- 
nal degrees of freedom - here projectively measuring the 
composition of the Raman dressed state - that allows 
the determination of p ca n- Figure 4 A shows the resulting 
time evolution of p can for several different initial A* all 
for Aj ^ 0; as expected, we observe oscillations of p can 
about q*A*j. Since At is small compared to our « 25 ms 
trap period, we expect the change of momentum to be 
dominated by Ap = —q*(A*f — A*), where the contribu- 
tion from the trapping force as A* changes is negligible. 
This translates into an oscillation amplitude Ap in both 
p can and m*v = p can — q* A^ of dressed atoms; the solid 
crosses in Fig. 3 show the amplitude of the sinusoidal os- 
cillations in p can as a function of Aj — A* w —^4*, proving 
that the synthetic electric field has imparted the expected 
momentum kick. 

We repeated the experiment with a nonzero q*A*f/h& 
0.35/cl, and observed, remarkably, the sinusoidal oscilla- 
tions in pcan were offset from zero as depicted in Fig. 4B. 
Pcan evolves smoothly about q* A^ ^ 0, illustrating the 
observed quantity is not the mechanical momentum mi>, 
which should oscillate about zero. We also measured mv, 
where v is the population-weighted average velocity of 
all spin components; while mv indeed oscillates about 
zero, the oscillation amplitude is smaller than that of 
Pcan- Given an increased effective mass m/m* w 0.4, the 
trap frequency v x along x was reduced by ^m/m* from 
that for undressed atoms, and the oscillation amplitude 
of tyiv should be reduced by m/m* from that of Pcan- 
Our results show a reduction factor of 0.38(4) in in 
agreement with the predicted m/m*, and the reduction 
factor of 0.30(2) in the momentum oscillation amplitude 
is slightly smaller than predicted. 

Conclusion We have demonstrated how a time- 
dependent, spatially uniform effective vector potential 
for low-energy dressed neutral atoms gives rise to a syn- 
thetic electric force. We have measured the effect of such 
forces for both trapped and untrapped atoms, and for 
the final state both dressed and undressed. While the 
electric forces generated here are spatially homogenous, 
this technique is generally applicable to create spatially 
varying forces. Indeed, since the effective vector poten- 
tial A* is parameterized by the Raman detuning 5, it 
can be locally patterned via suitable spatially inhomo- 
geneous magnetic fields or vector light shifts, giving rise 
to a nearly arbitrary range of time-dependent potentials. 
For azimuthal vector potentials, our method of measur- 
ing both the mechanical and canonical momentum from 
the induced electric fields can be used to identify the su- 
perfluid fraction of a BEC [T6] . or any other cold-atom 
systems. 
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funds from the DARPA OLE program, and the NSF 
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FIG. 4. Oscillating atoms in the trapping potential after application of a synthetic electric field pulse. (A)(B) Left panels: The 
vector potential is changed from q* A* = 0.75 hk\, (red circles), 0.25 hkz, (black circles), and q*A* = q*A*f « (green circles), 
respectively, all to q*A*f, and from q* A* = 0.75 hkh (blue circles) to q* Af = 0.35 hk^ (blue squares). The measured momentum 
for all circles is the canonical momentum p C an, and that for the square is the mechanical momentum mv. In (B), p can oscillates 
about q* A*f / while mv oscillates about zero. Right panels: Energy-momentum dispersion curves for uncoupled states (grey) 
and dressed states (colored). The arrows indicate oscillations of p can about q* Af for atoms in the lowest energy dressed state. 
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METHODS 

1 Dynamic change of effective vector potentials 

In our first measurement of synthetic electric fields, 
we observed the momentum imparted by the field kick, 
resulting from a change in the effective vector potential 



6 



from an arbitrarily chosen q* : A* to q* — ±2Hki, (see 
Fig. 3). In principle we could use any A*j and observe 
a momentum kick q*(A* — Af), however, in general the 
effective vector potential also depends on the strength 
of the Raman coupling Q,r. As a result, additional syn- 
thetic electric fields typically appear when Qr is adia- 
batically turned off. There are three specific cases for 
which A* does not depend on Qr: when the detuning 
5 = and 5 — >• ±oo. For the former case, only when 
< hkL is there no additional electric force dur- 
ing the removal of Qr, where the final atomic state is 
Imp = 0). For > hkL, the final atomic state is 

\rriF = ±1) with an additional momentum of ^■2hkp im- 
parted. Thus in our experiment we changed the vector 
potential from to q* A*j = ±2hk\ u by changing the 

detuning from Si to a large HSf = ^60El; the subsequent 
turnoff of then exerted no additional forces. 

2 Control of Raman detuning 

In all of our experiments, we set the Raman detuning S 
with small changes of the bias magnetic field away from 
resonance and hold the 2.32 MHz frequency difference 
between the Raman beams constant. Because all tem- 
poral changes in S lead to synthetic electric fields, bias 
magnetic field noise and relative laser frequency noise 
can lead to motion in the trap or heating. We phase 
locked the two Raman beams and observed no change in 
the heating, showing that relative laser frequency noise 
is not important in our experiment. However, our exper- 



iment is very sensitive to ambient magnetic field noise, 
here tied to the 60 Hz line. This noise gives rise to in- 
tractable dynamics of the canonical momentum of the 
dressed state, where S is held constant after the loading. 
We measured the field noise from the state-decomposition 
of a rf-dressed state (no Raman fields) nominally on res- 
onance and then feed-forward canceled the field noise. 
This reduced the ~ 0.2 jiT RMS magnetic field noise 
at 60 Hz by about a factor of 20, and remaining RMS 
field noise is ~ 0.03 fiT (including all frequency compo- 
nents up to « 5 kHz). All of our measurements were 
performed by locking to the 60 Hz line before loading 
into the dressed state. 



3 Momentum measurements of the dressed state 

The Raman-dressed state (no rf) is a superposition of 
spin and momentum components; its canonical momen- 
tum p can is the momentum of the spin mp = compo- 
nent. Experimentally, we fit the mp = density dis- 
tribution after TOF to a Thomas- Fermi profile [17] and 
identify p can as the center of the distribution. The me- 
chanical momentum of the dressed state mv was mea- 
sured by a population-weighted average over all three 
spin states including every pixel with discernable atoms 
in the image. This takes into account the modification 
of the TOF density distribution for all mp states due 
to interactions during TOF. Although interactions can 
exchange momentum between spin states, the total mo- 
mentum is conserved. 



